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Abstract 



This paper is concerned with highway traffic estimation using traffic sensing data, in 
a Lagrangian-based modeling framework. We consider the Lighthill-Whitham-Richards 



(LWR) model (Lighthill and Whitham 1955 Richards 1956) in Lagrangian-coordinates, 
and provide rigorous mathematical results regarding the equivalence of viscosity solutions 
to the Hamilton-Jacobi equations in Eulerian and Lagrangian coordinates. We derive 
closed-form solutions to the Lagrangian-based Hamilton-Jacobi equation using the Lax- 
Hopf formula ( Daganzo , 2005 Aubin et al. 2008 ) , and discuss issues of fusing traffic data 



of various types into the Lagrangian-based H-J equation. A numerical study of the Mobile 



Century field experiment ( Herrera et al. 2009 1 demonstrates the unique modeling features 



and insights provided by the Lagrangian-based approach. 



1 Introduction 



1.1 General background 

Highway traffic estimation and prediction are of pivotal importance for traffic operation and 
management. From an estimation perspective, it is desirable to have a substantial amount 
of information available. This has been made possible by emerging technology in fixed and 
mobile sensing and by the advancement of cyber-physical infrastructure, see Herrera et al. 



(2009) for a specific discussion. Depending on the characteristics of sensing devices and the 



underlying modeling framework, traffic sensing can be categorized as Eulerian sensing (such 
as using loop detector, virtual tripwire, and video camera) and Lagrangian sensing (such as 
using GPS and smart phone). 

This paper is concerned with using the Lighthill-Whitham-Richards model to estimate 
traffic states, subject to constraints imposed by sensing data that are large in quantity, and 
high in dimension. The LWR model is described by a partial differential equation (PDE), 
more specifically, by a scalar conservation law (SCL). For general mathematical background 
on SCL, the reader is referred to Bressan (2000). The LWR SCL is closely related to a 



Hamilton-Jacobi equation, via the Moskowitz function (Moskowitz, 1965), also known as the 
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Newell curve (Newell, 1993). It is shown in Evans (2010) that a weak entropy solution to the 
LWR SCL leads to the viscosity solution of the corresponding H-J equation; and vise versa. 



1.2 Viability theory and Lax-Hopf formula 

From a PDE perspective, additional information on the traffic stream provided by sensors 
imposes multiple initial/boundary and/or internal boundary conditions that may lead to non- 
existence of a well-defined solution, namely, entropy solution to the scalar conservation law; 
or viscosity solution to the H-J equation. Recently, Aubin et al. ( 2008| ) extended the viscosity 
solution of H-J equations to a solution class known as the Barron- Jensen/Frankowska solu- 
tions (Barron and Jensen, 1990 Frankowska, 1993). The generalized solution, called viability 



episolution (Aubin, 1991), is lower semi-continuous, and satisfies the initial/boundary and/or 
internal boundary conditions in an inequality sense. Such a relaxation, while excluding some 
nice mathematical properties of the solution, does allow the PDE to incorporate multiple 
value conditions and produce solutions that are well-informed of traffic states transmitted by 
the sensors. Based on this idea, a sequence of seminal papers, Claudel and Bayen (2010a|b 



2011) were created, which performed highway traffic estimation as well as data assimilation 



and inverse modeling. 



1.3 The LWR model in Lagrangian coordinate 

The LWR model is commonly formulated as a scalar conservation law in Eulerian coordinate 
(EC), i.e. with independent variables t (time) and x (location). A list of selected references 



Daganzo 


(2005 


2006 


); 


Friesz et al. 


(2012 


); 


Newell| ( 


1993) 



Another way of formulating the 
traffic dynamic is through Lagrangian coordinate (LC), which identifies t (time) and n (vehicle 
label). In contrast to EC, LC describes time evolution of traffic quantities associated with a 



moving vehicle. Such a framework was initially introduced by Courant and Friedrichs ( 1948 ) 
in the context of gas dynamics. LC was introduced to traffic flow modeling by |Daganzo ( 2006 ), 
and subsequently studied by Leclercq et al. (2007) and Yuan et al. ( 201 1[ ) . A more detailed 
review of the LWR model in both EC and LC is presented later in Section [2} 

The LWR PDE in LC can be derived from the Eulerian-based PDE; and vise versa. This 



procedure was given an intuitive illustration in Daganzo (2006), provided that all quantities 
of interest were continuously differentiable. In Leclercq et al. (2007), the authors related the 



equivalence of EC- and LC-based PDEs to an earlier work (Wagner, 1987) on gas dynamics. 



However, we point out that Wagner (1987) studied an Euler equation, which is a system of 
conservation laws; results regarding such a system do not apply immediately to the scalar 
conservation law of interest. In this paper, we provide a rigorous proof of the equivalence 
between the two systems, in the context of Hamilton-Jacobi equation and viscosity solutions. 
The result will apply to the LWR conservation laws even if the solutions are not continuous, 



let alone differentiable. Technical result is presented in Section 2.3 

As mentioned before, solution to the H-J equation has been extended to treat multiple 
value conditions, using viability theory (Aubin, 1991). In this paper, we apply the frame- 
work of viability theory to the Lagrangian-based Hamilton-Jacobi equations. The goal is 
to explore the unique modeling and computational advantages of the LC in highway traffic 
estimation. In particular, using the Lax-Hopf formula, we derive explicit solutions to the 
LC-based H-J equation with various value conditions including initial/intermediate condition, 
upstream/downstream boundary condition, internal boundary condition, and a combination 
of the above. This is presented in Section |4j As pointed out by recent studies Leclercq et al. 
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(2007) and Yuan et al. (2011), LC-based conservation law has certain computational advan- 



tage over EC-based ones, in the context of Godunov scheme (Godunov, 1959). Namely, the 
Godunov scheme reduces to an upwind scheme for LC-based conservation law, since the cor- 
responding fundamental diagram is monotonically increasing. As we show in this paper, when 
solved with the Lax-Hopf formula, the LC-based H- J equation also reduces the computational 
complexity, compared with the EC-based H-J equation. Finally, we test the framework of 
Lagrangian-based traffic sensing and estimation with dataset obtained from the Mobile Cen- 



tury field experiment ( Herrera et al. , 2009 ) . The Lagrangian-based approach provides traffic 
information associated with each moving vehicle, such as trajectory and speed. Such informa- 
tion are relatively difficult to obtain from the Eulerian-based approach: one has to perform 
partial inversion of the Moskowitz function. 

1.4 Organization 

The rest of this paper is organized as follows. In Section [2j we briefly review the Lighthill- 
Whitham- Richards model in both Eulerian and Lagrangian coordinates. A rigorous equiva- 



lence result regarding the EC-based and LC-based H-J equations is presented in Section |2.3 
Section [3] introduces the notion of viability episolution and the generalized Lax-Hopf formula. 
Data of different sources are discussed and integrated into the Lagrangian PDE. In section [4j 
we derive explicit solutions to the Lagrangian-based H-J equation, using the Lax-Hopf formula. 
In Section [5j we apply the Lagrangian-based computational paradigm to estimate real- world 
highway traffic stream. 

2 The LWR model in transformed coordinates 

In this section, we review the Lighthill-Whitham-Richards model in both Eulerian and La- 
grangian coordinates. The transformation between these two coordinate systems is made 



through a function inversion under minor assumption. In section 2.3, rigorous result on the 
equivalence of the two systems is provided, using Hamilton-Jacobi equation and the notion of 
viscosity solution. 

2.1 The LWR model in Eulerian coordinates 

The classical LWR PDE reads 

Pt{t, x) + f x (p(t, x)) = (t, x) E [0, +oo) x [0, L] (2.1) 

where the subscripts denote partial derivative. The scalar conservation law above describes 
the temporal-spatial evolution of average vehicle density p(t, x) and average vehicle flow 
f(p(t, x)). The density-flow relation is articulated by the fundamental diagram /(•) 

f(p) = pv(p) pE [0, Pmax ] (2.2) 

where p m ax is jam density. The average vehicle speed v(p) E [0, v max ] is a decreasing function 
of density; v max denotes the free flow speed. The fundamental diagram is assumed to be 
concave with maximum M attained at p* , where M is the flow capacity. 

We introduce the Moskowitz function N(- , •), also know as the Newell-curve, such that 

N t (jt,x) = f[p(t,x)), N x (t,x) = -p{t,x) a.e. (2.3) 
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Hereafter, "a.e." stands for "almost everywhere". N(t, x) measures the cumulative number of 
vehicles that have passed location x by time t. It has been shown in a number of circumstances 



that if pit, x) is the weak entropy solution to (2.1 ), then the corresponding Moskowitz function 



is the viscosity solution to the following Hamilton-Jacobi equation 



N t (t, x)-f(- N x (t, x)) 







(2.4) 



Note that viscosity solution to ( |2.4[ ) must be Lipschitz continuous, but not necessarily 
continuously differentiable due to the presence of shock waves in p(t, x). There exists, however, 



more general solution classes such as the viability episolution (Aubin et al. 2008), which is 
lower semi-continuous. 



2.2 The LWR model in Lagrangian coordinates 

In Lagrangian coordinate system, a free variable n is used to identify a particular vehicle. Note 
that in a continuum model, n is treated as a real number. The coordinate transformation from 
(t, x) to (t, n) is made by inverting the following Moskowitz function 

n = N(t, x) (2.5) 

Throughout this article, we assume that the vehicle density is uniformly positive, i.e. there 
exists 5 > such that 



p(t, x) > 5 for all (t, x) G [0, +oo) x [0, L] 



(2.6) 



This assumption is not restrictive given the argument that if the density vanishes in certain 
road segments, then the domain of study can be separated into several subdomains, with each 
one satisfying (2.6). 



For each t fixed, (2.6) implies that N(t, •) is a strictly decreasing function of x, whose 

x = X(t, n) (2.7) 



inverse is denoted by X(t, •). We have 



where X(t, n) represents the location of vehicle labeled n at time t. Identities (2.5) and (2.7) 
define the coordinate transformation. We have the following 



Xt(t, n) = v(t, n), X n (t, n) 



s(t, n) 



a.e. 



21 



where v(t, n) and s(t, n) denote the average speed and spacing around vehicle n at time t, 
respectively. Introducing the spacing-velocity relationship 



11,(8) = f 



S£[0,V max ] S G [1/ Pmax, 1/ $] 



we present the Hamilton-Jacobi equation in Lagrangian coordinate 

X t (t,n)-^(-X n (t,n)) = 



(2.9) 



(2.10) 
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2.3 Equivalence between (|2.4[) and (|2.10[) 



In this section, we provide rigorous mathematical analysis on the relationship between the 
two H-J equations in EC and in LC. Our strategy is to show that if N(t, x) is the viscosity 



solution to (2.4), then its partial inverse as defined by (2.5) and (2.7), is the viscosity to 



(2.10). The converse will hold similarly. We begin with the definition of viscosity solution to 



Hamilton-Jacobi equation of the form 



u t + H(Vu) = 



(2.11) 



where the unknown u(t, x) G M. m ; Vit is the gradient of u with respect to x. In what follows, 
C, C 1 denotes the set of continuous and continuously differentiable functions, respectively. 



Definition 2.1. A function u G C(Q) is a viscosity subsolution of (2.11) if, for every C 1 



function ip = ip(t, x) such that u — cp has a local maximum at (t, x), there holds 

tp t (t, x) + H(Vip) < (2.12) 



Similarly, u G C(fi) is a viscosity supersolution of (2.11) if, for every C 1 function ip = ip(t, x 



such that u — ip has a local minimum at (i, x), there holds 

<pt(t, x) + H(V<p) > 



(2.13) 



We say that u is a viscosity solution of (2.11) if it is both a supersolution and a subsolution 
in the viscosity sense. 



Remark 2.2. If u is a C 1 (r2) function and satisfies (2.11) at every x G then u is also a 
solution in the viscosity sense. Conversely, if u is a viscosity solution, then the equality must 
hold at every point x where u is differentiable. In particular, if u is Lipschitz continuous, then 



it is almost everywhere differentiable, hence (2.11) holds almost everywhere in Q. 



The next theorem establishes the equivalence between viscosity solutions to (2.14) and 



(2.15), where /(•) and ifj(-) satisfy (2.9). 



N t (t, x )-f( K -N x (t, x)) 



(2.14) 

X t (t, n)-il>(-X n (t, n)J = (2.15) 

Theorem 2.3. Assume that N(t, x), (t, x) G C (— oo, +oo) x W 1 , is a viscosity solution 
to (2.14), furthermore, assume that the density is uniformly positive, i.e. p(t, x) > 5 > 
0, for all (t, x) G $7. Then function X(t, •) obtained by inverting N(t, •) is a viscosity solution 



to $2.15] ). 

Proof. By assumption, N(t, •) is strictly decreasing with 

5\xi-X 2 \ < \N(t, Xi) - N(t, X 2 )\ < Pmax\x\ 

then X(t, ■) is also strictly decreasing with 

l/Pmax\ni ~ 712 1 < \X(t, Til) - X(t, 71 2 ) | < l/5\ni-n 2 



X2\ 



for all xi, X2 



for all ni, ri2 



We start by showing that X(-, •) is a subsolution. Indeed, given any C 1 function Y = 
such that X — Y has a local maximum at (to, no), without loss of generality, we 
X(to, no) — Y(to, no) = 0, take the plane t = to (see Figure [T]). 



(2.16) 

Y(t, x) 
assume 
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Figure 1: Graphs of X(to, •) and Y(to, •) 



Since X — Y attains a local maximum at (to, no), by (2.16) there must hold Y n (to, no) < 0. By 
continuity, there exists a neighborhood S7i of (to, no) such that Y n (t, n) < 0, for all (t, n) £ 
ill. Then, we may define M(t, •) to be the inverse of Y(t, •) in Q±. Obviously, M(t, n) E 
C 1 (Oi), and N — M attains a local maximum at (to, X(to, no)). Using the fact that N(t, x) 
is a viscosity solution and applying (2.12), we deduce 



M t (t,x) < f(-M x (t,x)^J 



(2.17) 



Differentiating with respect to t the identity Y(t, M(t, x)) = x, and using (2.17), we have 



Y t + Y n M t > Y t + Y n f 



Y t + Y n f{ 



(-~) = Y t -i;(-Y n ) (2.18) 



In the above deduction, we have used the fact that M(t, •) and Y(t, •) are both C 1 , and are 
inverse of each other. Therefore, differentiating n = M(t, Y(t, n)) w.r.t. n yields 1 = M x • Y n . 

Since Y is arbitrary, (2.18) implies that X(t, n) is a subsolution. The case for supersolution 
is completely similar. □ 

Remark 2.4. Similar proof can be used to show the converse: given a viscosity solution X(-, •) 
to (2.15), N(t, •) obtained by inverting X(t, •) is a viscosity solution to (2.14). 



3 Numerical algorithm and value conditions 

In this section, we focus on the numerical scheme for the Hamilton-Jacobi equation 

d t X(t,n)-7p(-d n X(t,n)) = (3.19) 
in the presence of initial, boundary and internal boundary conditions. In order to avoid the 



issue of non-existence of solution, we adapt the notion of viability episolution (Aubin et al. 



2008 



2008 



, and proceed with a variational method known as the Lax-Hopf formula (Aubin et al. 
Claudel and Bayen[ [2010^ . 



3.1 Viability episolution to the Hamilton-Jacobi equation (3.19) 



We consider a domain [0, T] x [Ni, N2] for equation (3.19), where T > and N2 > N\ > 0. 



In the following definition, we make precise what we mean by value conditions, concerning 
d3l9b. 
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Definition 3.1. A value condition C(-, •) : O C [0, T]x[Aq, iV 2 ] 



is a lower- semicontinuous 



One may extend a value condition C(-, •) to the entire domain by assigning C(i, n) = +00 
whenever (t, n) £ O. Such a convention allows one to compare and manipulate value conditions 
with different domains. In order to articulate the Lax-Hopf formula, we introduce the concave 
transformation ip*{-) °f the Hamiltonian ip(-) 

4>*(p) = SU P w( s ) — ps\ 

se[l/p m ax,+Oo) 

The following Lax-Hopf formula provides semi-analytical representation of viability solution, 
given a value condition. 



Theorem 3.2. The viability episolution to (3, lty associated with value condition C(-, •) is 
given by 

X c (t,n)= inf (C{t-T, x + Tu) + Ti;*(u)) (3.20) 

(u,T)eflom(iii*)xi + 



□ 

We indicate, by a subscript, dependence of the viability solution on its value condition. So 



Proof. See Aubin et al. (2008) 



that the solution to (3.19) reads Xc(t, n). As a consequence of formula (3.20) the following 
property of inf-morphism holds. 

Proposition 3.3. (inf-morphism property) Let C(-, •) be the point wise minimum of 
finitely many value conditions, 



C(t, n) 



min d(t, n) for all (t, n) G [0, T] x [Aq, N 2 ] 

i=l,...,m 



Then 



X c {t, n) 



min X Ci (t, n) 
i=l,...,m 



(3.21) 



The inf-morphism property enables the H-J equation to incorporate one or more value 
conditions; it can also decompose a complex problem involving multiple value conditions into 
smaller subproblems, each one with a single condition. 



3.2 Value conditions for continuous solutions 

The value conditions described in the previous section can be viewed as the mathematical ab- 



straction of real-world measurements, obtained from Eulerian and Lagrangian sensing. Claudel 



and Bayen ( 2010a|b ) have utilized the H-J equation in Eulerian coordinates to incorporate 
sensing data. In the rest of this paper, we not only extend their framework to Lagrangian 
coordinates, but also establish sufficient condition that guarantees the equivalence between 
these two systems, in the presence of a single value condition. 



Consider again the two H-J equations (2.14) and (2.15) in EC and LC. We make two 



assumptions on the value condition C, as follows 

(Al) The domain of C is a continuous curve parametrized by r G [r m j n , r max ] 

Dom(C) C [0, T] x [N u N 2 ] (respectively [0, T) x [X 1} X 2 ]) 
Dom(C) = (t(r), n(r)) (respectively (t(r), x(r) 
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(A2) C(t(-), n(-)) is a continuous function on [r m j n , T max ]. 

Given an arbitrary continuous value condition C(-, •) in the Eulerian domain (e.g. measure- 
ment of a loop detector), we need to determine how to fuse such datum into the Lagrangian 
based PDE. The next proposition provides an answer. 

Proposition 3.4. (Sufficient condition for equivalence of value conditions) 

Let C E (t(r), x(t)) and C L (t(r), n(r)), r S [r m j n , T max ] be two value conditions for (2.14) 



and (2.15), respectively. Then the solutions to (2.14) and (2.15) satisfying each value condition 



in the equality sense are equivalent if 



n(r) = C E {t(r),x(r)), x(r) = C L '(i(r) , n(r)) 



Proof. Let N^eU, x) be the solution to ( |2.14[ ) satisfying condition C E , let X C h(t, n) be the 
solution to (2.15) satisfying condition C L . For eac h t, s ince N C E(t, •) is strictly increasing, we 
denote its inverse by N~ E (t, •). Then by Theorem 



equation (2.15). 



2.3 



iV" 



is a valid solution to the HJ 



On the other hand, for every r £ [r m j n , r„ 



iV-i(i(r), n(r)) = N^(t(r), C E (t(r), x(r))) = (t(r), AT cE (t(r), x(r))) = s(r) 

(3.22) 



(3.22) implies that N CE (-, ■) satisfies value condition C L (t(r), n(r)) and thus is the unique 
solution to (2.15) associated with value condition C L . We conclude N^ E (t, n) = X C L(t, n). 

□ 



Proposition 3.4 suggests that, by simply invert the Eulerian (Lagrangian) value condition 



C(t(r), •), one can obtain its dual version in LC (EC). The dual value condition yields solution 
that is equivalent to the original solution. 



4 Explicit solution with piecewise affine value conditions 



In this section we apply the Lax-Hopf formula (3.20) to for an explicit solution in the La- 
grangian coordinates. To that end, we assume that the value conditions are piecewise affine 
(PWA). This assumption is justified on the ground of (1) an easily constructible explicit so- 
lution representation, with PWA conditions; and (2) standard result on PWA approximation 
of any piecewise continuous functions or curves. 



4.1 Piecewise affine value conditions 

We start with articulating piecewise affine value conditions, which includes initial/intermediate, 
upstream, downstream and internal conditions. 

Definition 4.1. (PWA initial/intermediate condition) 

Set t = to > 0, given real numbers Sj > 0, rtj, i G {1, . . . ,rrii n i}, the j th affine component 
of initial/intermediate condition is 

Cini(t, n) = -Sjn + dj, (t, n) £ {t } x [nj, n j+1 ] (4.23) 

To ensure continuity, we require 

i-i 

d j = s 3 n 3 ~ 5Z( nm ~ n ')' = 2 ' • • ' ' mini 

l=i 



S 



Definition 4.2. (PWA upstream boundary condition) 

Fix n = N\, given real numbers > 0, t\ i G {1, . . . , m up }, the j th affine component of 
upstream boundary condition is defined as 

Ci p {t, n) = vH + b>, (t, n) £ [t j , t J+1 ] x {N^ (4.24) 

To ensure continuity of upstream boundary condition, we set 

3-1 

V = - v Hi + Y,(t l+1 j = 2,...,m up 

i=i 

Definition 4.3. (PWA downstream boundary condition) 

Fix n = N2, given real numbers V{ > 0, ti, i G {1, . . . , rridown}, the j th affine component of 
downstream boundary condition is defined as 

C L™(*> n ) = V + b J> (*» n ) G x {^2} (4.25) 

j'-i 

b j = ~ U i h + y^fe+l - *0 U J> i = !> • • • > "^oton 

z=l 

Definition 4.4. (Affine internal boundary condition) 

Given real numbers a, (3,t m i n , t max , n m i n , n max , r > 0, i/ie affine internal boundary con- 
dition is defined as 

C%nt\b} Tl) — j3 -\- Oi {t t m in) t G \tmim t max \i Tl — f^min ~\~ r(t ^min) (4.26) 

Recall that the domain of our consideration is [0, T] x [N\ , N2] , thus the upstream/downstream 
boundary conditions refer to (part of) the trajectories of the first and last car within our scope. 

4.2 Explicit formulae for viability episolutions 



In the presence of piecewise affine (PWA) data, the solution of Lagrangian equation (2.15) 
with a piecewise affine Hamiltonian can be computed explicitly, as given by the Lax-Hopf 
formula (3.20). We denote s m i n = l/p max , s* = 1/p*, let k > 0, define 

Ms) = (*'' ^[w, S *] (427) 

[Vmax, S G [S , +Oo) 

Notice that this fundamental diagram corresponds to a triangular density-flow relationship. 
Moreover, the concave transformation of ip(-) reads 

^*(u) = s*(k-u) ue[0,k] (4.28) 



Proposition 4.5. Given each affine value condition defined in (4-23)-(4-26), and Hamiltonian 
(4-21), the solution to the Lagrangian Hamilton- J acobi equation (2.15) are respectively 

1. Initial/intermediate value problem 
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ifsj < s*, 



—Sj n + dj + (t — to)k Sj, 

rij + k(t — to) < n < rij + i + k(t — to); 

—Sj rij + dj + (t — to) s* k — s* (n — rij), 
< n — rij < k(t — to). 



(4.29) 



if Sj > s*, 



—Sj nj + i + dj + (t — to) s* k — s* (n — rij+i), 
< n — rij + i < k(t — to); 

-Sj n + dj + v max (t - to), 
rij < n < 



(4.30) 



Upstream boundary value problem 



V#' +1 + 6? + s*(fc(* - t j+1 ) - (n - JVi)), 
< n-Ni < k(t-V +1 ); 



X J up (t, n) = < 



v j t + V - (n-AM— , 
k 

^max{0, k(t - t j+1 )} <n-N x <k{t- P). 



(4.31) 



3. Downstream boundary value problem 



X down(^ X ) = V 3 t 3+ l + h 3 + (* ~ (*> x ) G fe+l> +°°) X {^2} 

^. Internal boundary value problem 

n , / W Tlmin kt -\- V t m i n 

P + a[ : t m 

\ r — k 

H 17"min ^ f \t tmin), 
k(t t max ) "C. Tl fl max , 
Tl 1T"min — k(t imm)j 



X int (t, n) 



P H a + S kit- t m i n ) 

^ ti fi m i n <C r(t t m { n ) and n ^ Ti max , 



ft ~\~ OL (tmax train) \t t max ) S k S (jl Tlmax), 

k ^ ti fi max ^ kit t max ), 



(4.32) 



(4.33) 



Remark 4.6. (4-33) is well defined even if k = r and r = 0. 
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The proof of Proposition 4.5 is straightforward. A simple verification reveals that the 
formulae for upstream, downstream and internal boundary conditions (4.31 )-(4.33) can be 
merged into one, namely (4.33). In other words, the upstream and downstream conditions 
can be treated as spacial cases of internal boundary conditions. This is because the wave 
propagation speed is always non-negative, and any value condition can only affect the solution 
at region with a larger n. This coincides with the intuition that a car cannot be affected by 
others behind it. 

In order to compute the PDE solution with more complicated value conditions, for example, 



the one including multiple data, we invoke the inf-morphism property from Proposition 3.3 
and take the point wise minimum over all solutions obtained from (4.29), (4.30) and (4.33). 

The Lax-Hopf formula is a grid-free numerical scheme in the sense that the computational 
procedure does not rely on a two-dimensional grid, as opposed to finite-difference schemes. 
In addition, the solution algorithm is highly parallelizable: the problem of fusing multiple 
value conditions into one PDE can be decomposed into independent sub-problems. In each 
sub-problem, the solution can be solved efficiently with formulae (4.29)-(4.33). 



5 Numerical Study 



In this section, we conduct a numerical case study using traffic data obtained from the Mobile 



Century field experiment (Herrera et al. 2009). The dataset includes vehicle trajectory data 



and cumulative passing vehicles measured at different locations along a highway segment. We 
utilize a subset of the Mobile Century dataset and fuse them into the Lagrangian Hamilton- 
Jacobi equation. The goal is to perform the estimation of Lagrangian quantities such as vehicle 
trajectories and velocities. 



5.1 The Mobile Century field experiment 

On February 8, 2008, an experiment in the area of traffic monitoring was launched between 
9:30 am to 6:30 pm on freeway 1-880 near Union City in the San Francisco Bay Area, California. 
This experiment involved 100 vehicles carrying GPS-enabled Nokia N95 phones. The probe 
vehicles repeatedly drove loops of 6 to 10 miles in length continuously for 8 hours. 

Carried by each probe vehicle, each on-board smart phone stored its position and velocity 
every 3 to 4 seconds. In addition to the cell phone data, the experiment also collected data 
on cumulative passing vehicles at different locations along the highway. The reader is referred 
to 



Herrera et al. (2009) for more details on experimental design and data description. 



5.2 Details of the numerical study 

The freeway segment of our study is a 3.45 mile stretch of 1-880 North Bound, between 
postmile 23.36, and postmile 26.82. We utilize two types of data: (i) the cumulative vehicle 
count, obtained via loop detector station 400536 (postmile 23.36), which counts the number of 
passing vehicles every 30 seconds; (ii) vehicle trajectories, recorded by on-board smart phones 
every 3 to 4 seconds. Our study spanned one hour, from 11:30 am to 12:30 pm, and involved 
a total traffic volume of approximately 5000 vehicles. We utilized data sent from 97 probe 
vehicles, The vehicle trajectories are plotted in Figure [2] In order to identify the 'label' of 
those probe vehicles, we processed the data on the cumulative vehicle count obtained from 
Station 400536. 
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We assume that the fundamental diagram tp(-) in (2.15) is piecewise affine. Such a fun- 
damental diagram is calibrated with data collected form the same experiment (detail of the 
model calibration is omitted for brevity). The key parameters are estimated as follows. 

Vmax = 31.5 (meter /second), p m ax = 0.50 (vehicle /meter), p* = 0.055 (vehicle /meter) 

The resulting fundamental diagram reads 



Ms) 



5.3 Numerical results 



1.95 s if 2.00 < s < 18.15 

(5.34) 

31.5 if s > 18.15 



We reconstructed the traffic state between postmile 23.36 and 26.82 for a period of one hour, 
using the numerical method discussed in Section |4| Among the 97 vehicle trajectories, we 



used 44 of them as training data. In other words, the viability episolution to equation (2.15) 
was computed with the Lax-Hopf formula and 44 internal boundary conditions. Notice that 
the penetration rate of probe vehicles is around 0.88%, given the total traffic volume in this 
one-hour period. Recall that the on-board smart phone recorded vehicle location every 3 to 4 
seconds. However, in our actual computation, we sampled the vehicle location at a much lower 
frequency, namely at every 30 to 90 seconds. As we show later in this section, the numerical 
performance of the model remains robust under such crude measurements. 




Figure 2: Recorded trajectories of 97 probe vehicles. 



Figure [3] displays the solution of the H-J equation (3.19) based on 8 vehicle trajectories 
(internal boundary conditions). It should be noted that the viability episolution satisfies the 

that is 



value conditions only in an inequality sense (Aubin et al. 

X c (t, n) < C(t, n), 



2008), 



(t, n) G Dom(C) 
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If the strict inequality holds, the value conditions and the model are said to be incompati- 
ble. The incompatibility is due to either error in measurements or inaccuracy of the model 
itself. A class of problems derived from incompatibility, namely data assimilation and data 
reconciliation, were discussed in Claudel and Bayen (2011). 




'11 :33:20 7950 Vehicle Label 
Local Time 



Figure 3: Solution to the Hamilton- Jacobi equation (3.19) with 8 (internal) boundary condi- 
tions. The solid lines are plots of vehicle trajectories. 

The obtained solution X(t, n) describes the time evolution of the location of vehicle labeled 
n. In other words, it can be used to reconstruct vehicle trajectories. We fixed several values 
of n, and plotted the curves of X(-, n). Two examples are shown in Figure [4] and [5| For 
comparison purpose, we also included the actual vehicle trajectory reported by the smart 
phone in the same figure. 

Recall that 

d t X(t, n) = i/>(s(t, nj) = v(t, n) (5.35) 

By numerically differentiating X(t, n) with respect to t, one obtains the velocity field v (i, n). 
Figure [6] shows the velocity field involving vehicle labeled 7650 to 7950. This result was 
obtained by differentiating the Moskowitz function displayed in Figure [3} From Figure |6j we 
clearly observe a time- varying congestion level experienced by the drivers: vehicles experience 
first free flow traffic, then congested traffic, and finally less congested traffic. This observation 
is confirmed in Figure [2j where a similar driving condition is shown for all probe vehicles. 

Finally, we performed the task of travel time estimation using the Lagrangian approach. 
The travel times were estimated from the vehicle trajectories, which were again obtained by 
the viability episolution X(t, n). In order to examine the accuracy of our estimation, we 
compared the actual travel times of the 97 probe vehicles with the aforementioned estimated 
travel times. The result is summarized in Figure [7j 
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Figure 4: Reconstructed and measured trajecto- Figure 5: Reconstructed and measured trajecto- 
ries of vehicle #8685 ries of vehicle #11266 
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Figure 6: Vehicle velocity field (in meter/second) obtained from (5.35) 



6 Conclusion 

This paper applies the Lagrangian-based first-order hydrodynamic traffic model to perform 
highway traffic estimation, utilizing both Eulerian and Lagrangian sensing data. We briefly 
discuss the Lighithill-Whitham-Richards model in both Eulerian and Lagrangian coordinate 
systems, and provide rigorous mathematical results on the equivalence of Hamilton-Jacobi 
equations in these two coordinate systems. We derive closed-form solutions to the LC-based 
Hamilton-Jacobi equation, and discuss the issue of fusing traffic data of different sources into 
the H-J equation via the notion of viability episolution. A case study of the Mobile Century 
project demonstrates the unique modeling features and insights offered by the Lagrangian- 
based PDE. 

Despite several similarities between the EC-based and LC-based approaches such as struc- 
ture of the PDE and solution method, the LC has the advantage of providing information 
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Error (second) 
Figure 7: Error of travel time estimation 

associated with a given vehicle. Such a feature needs to be further explored and engineered 
to facilitate the development of cyber-physical infrastructure, including the mobile internet. 
In addition, extension of LC-based PDE to a network setting is also an important aspect of 
future research. 
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